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Theorems are presented concerning conditions for the input-output 
stability of interconnected dynamical systems. Results in the area of 
input-output stability are often partitioned into two categories: 
small-gain type-results and passivity-type results. The main theorem 
given here does not fall into either of these categories, but is most closely 
related to the passivity-type results. The theorem involves a new class 
of interconnection operators that is a substantial generalization of the 
familiar set of nonnegative operators defined on a space of vector- 
valued functions. 

I. INTRODUCTION 

In this paper, theorems are presented concerning the input-output 
stability of interconnected systems. 1 " Results in the area of input-output 
stability are often partitioned into two categories: small-gain type results 
and passivity-type results. The main theorem given here does not fall 
into either of these categories, but is most closely related to the passiv- 
ity-type results. The theorem involves a new class of interconnection 
operators that is a substantial generalization of the familiar set of non- 
negative operators defined on a space of vector-valued functions. 

The mathematical model considered throughout the paper is described 
in Section II, and results of a general nature concerning the model are 
given in Section III. The case in which the interconnection operator has 
a certain matrix representation is treated in considerable detail in Sec- 
tion IV. In Section 4.5, a specific example is given of a stable intercon- 
nected system for which the interconnection matrix does not meet the 
nonnegative-definiteness requirement of the criterion given in Ref. 7, 
which contains the most pertinent earlier stability result. 



* This paper was presented at the 1978 IEEE Symposium on Circuits and Systems (New 
York, May 17-19, 1978). 

t For background material in book form concerning intput-output stability, see, for ex- 
ample, Refs. 1-4. Interconnected systems (which are systems whose natural or artificial 
decomposition into subsystems plays a prominent role in their mathematical analysis) 
have been considered by many researchers. See, for example, Refs. 3, 5, 6, and 7. Although 
some interesting and significant results have been obtained concerning the stability of 
interconnected systems, the theory is very much in its initial stages of development. 
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The main purpose of the paper is to introduce a new concept that is 
believed to be useful. No attempt is made to present the sharpest possible 
stability results that the concept can be used to obtain. 

II. THE MODEL 

2.1 Preliminaries 

Let K denote a real linear space that contains a normed linear inner- 
product space L with inner product (•,•) and norm | • | related by | / 1 = 
(/»/) ^ 2 f° r f €zL. (Of particular interest to us is the case in which L is the 
set Li of all real Lebesgue square-integrable functions defined on the 
half line [0,°°) with the usual inner product, and K is the "extended" set 
E2 of real functions defined on [0,°°) such that each function is square 
integrable on [0,t] for any nonnegative number t.) 

For each r > 0, let P T denote a linear mapping of K into L (e.g., if K 
= E 2 , let P T be defined by (P T f)(t) = f(t) for t £ [0,t] and f(t) = for t 
> t, where / is an arbitrary element of E2). 

Let K, L, and P T be such that (0 g £ L if and only if g £ K and sup T 
|P T ^| < », (») \g\ = sup T \P T g\ forg £ L, and (Hi) (P T f,g) = {P T f,P T g) 
and \P T f\ < l/l for / and g in L and r > 0. 

We let L n and K n , in which n is any positive integer, denote the n-fold 
Cartesian product of L and K, respectively. The norm of an element h 
— (h\,h,2, ■ . ■ ,h n ) of L" is denoted by \h\ and is defined by \h\ = (2 t - 
|/i,| 2 ) 1/2 . 

It is assumed that L contains n elements e\,e2, ■ ■ ■ ,e n such that |e; | 
= 1 for each i and (e,,e ; ) = for i ^ j* 

We say that an operator T that maps K into itself (i.e., an operator 
T in K) is causal if and only if P T T = P T TP T on K for all r > 0. 

2.2 The basic equations 

Throughout the papei, attention is focused on an interconnected 
system governed by 

n 

Xi + E AijBjXy =y t , i = 1,2, . . . ,n, (1) 

>-i 
in which (A.l): xi and y^ belong to K for all i, and Ay and 5; are causal 
operators in K for all i and /. 

In (1), each Bj is associated (sometimes somewhat indirectly) with a 
subsystem, and the Ay ordinarily take into account the way in which the 
subsystems interact. Typically, it is not difficult to show the existence 
of a solution Xi, X2, . . . ,x n of (1) for any given y it y2, ■ . . ,y n under some 
weak additional hypotheses. (Successive-approximation type arguments 
of the kind commonly used in connection with nonlinear Volterra inte- 
gral equations often suffice.) 



* This assumption is used only in Section IV. 
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We assume that {A.2): each Bj is nonnegative in L, in the sense that 
each Bj maps L into L and there exists a nonnegative constant a such 
that 

{BjW,w)>a\w\ 2 (2) 

for w E L and all ;'. It is assumed also that (A.3): each Ay maps L into 
itself, and there is a positive constant y such that 

\Aijw\<y\w\ (3) 

for w E £ and all i and /.* 

It is often convenient to write (1) in the form 

x + ABx = y, (4) 

in which x = (xi,x 2 , • • • ,x n ), V = (yiJ* ■ ■ ■ Jn), and A and B are the 
mappings of K n into K n defined by (A/), = 2/Ay£ and (B/), = B,-/,- for 
all/ E # n and each i. 

2.3 Definition of stability 

We say that (4) is L-stable if and only if y E £ n implies that x E L n 
with |x | < p(\y |) for some nonnegative continuous function p that de- 
pends only on A and B and is defined on the nonnegative reals such that 
P(0) = 0. 

III. S f , AND THE MAIN THEOREM 

In the following definition and hypothesis, is a nonnegative num- 
ber. 

Definition of S : Sp is the set of operators H in L n with the following 
property: For each v E L n such that |u| ^0, there is an index k such 
that | u* | ^0and(t;fc,(//D)/e)>/3M 2 . 

The definition of S is related to one of the equivalent definitions of 
a "P-matrix." 9 
H.J: If (I ■ 0, there is a positive constant 5 such that 

|J3;iu| < d|w| (5) 

for w E ^ and all /. 

Our main result is the following. 
Theorem 1: Let H.l (and A.l through A.3) hold. Then (4) is L-stable 
ifAESp with a + P>0. 



* In order to focus attention only on essentials, we are proceeding with some assumptions 
that are stronger than necessary. It will become clear that (2) and (3) (and also (5) of Section 
III) could have been replaced with somewhat weaker inequalities (similar, for example, 
to some of those used in Section 5.3 of Ref. 8). Similarly, if for example there are positive 
constants a ; such that (BjW,w) > a; | w \ 2 for w 6 L and all ;", and if A in (4) is represented 
by an nXn matrix in the sense of Section IV with [I + a diag (a,)] invertible, then it is clear 
that x satisfies an equation similar to (4) in which each BiXj is replaced with {fijxi - afitj) 
and A and y are modified accordingly. Consideration of such a modified equation often 
enables a useful trade-off to be made between requirements on A and the degree of posi- 
tiveness of the Bj. 
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Proof of Theorem 1 : Let H.l and A. 1 through A.3 be satisfied, let y G 
L n , let x £ K" be a solution of (4), and assume that A £ S^ with 
a + (8 > 0. 

Suppose first that (3 = 0, in which case a > 0. Let t > be arbitrary. 
Using (1) and the casuality of the Ay, we have 

PrXi +P T Z AijPrBjXj = P r V ( (6) 

i 
and 



(P&XifrXi) + (P T B iXi ,P T Z AijPrBjXj^ = (PrBiX^Pryd 



(7) 



for i = 1,2, ... ,n. Since A G S^, and {.P T B i Xi,P r l l jAijP T BjXj) = 
(PjBiXi^jAijPrBjXj) for each i, there is an index k such that 
(P r BkXk,P T ^jAkjP T BjXj) > and hence such that 

(P T B k x k ,P T x k ) < (P T B k x k ,P T y k ). (8) 

By the Schwarz inequality and the fact that {P T BkXk,P T yk) = 
(P T B k P T x k ,P T y k ) = (BkPrXk.Pryk), (P T B k x k ,P T y k ) < \B k PrX k \ ■ \P T y k \. 
Therefore, using (P T BkXk,P T Xk) = (BkPrXkfrXk) as well as (2), (5) and 
(8), we have 

a|P r x,|2<5|P T ;c,|.|P T ;>/ fe | (9) 

and consequently, with c = 5a -1 , 

\P T x k \ <c\P T y k \. 

The argument given above shows that \PjXk \ ^ c \y \ for some k (which 
might depend on x and t). Let J denote any nonempty proper subset 
of {1,2, . . . ,n\ with the following property. For j £ J, there is a constant 
Cj that depends only on c, 5, and y such that |P T *y| — c y -|y|. Using 
(1), 

x t + E AijBjXj = yi- E AijBjXj, i € J. (10) 

The left side of (10) is basically the same in form as the left side of (1). 
With r the number of elements contained in J, let the elements of 
({1,2, . . . ,n\ - J) bej h J2, ■ ■ ■ ,j( n -r) ordered so that;'i <j 2 < • • • <j( n -r)- 
With respect to that ordering, let the mapping of K (n ~ r) into itself as- 
sociated with (10) that corresponds to A be denoted by Aj. Since A.3 
holds, each A, ; maps the zero element of L into itself, and it is a simple 
matter to verify that Aj belongs to, so to speak, Sp with n replaced with 
(n — r). Thus, by the argument given above, we find that there is an index 
/ € J such that 

|P T x,| < c|p T (y, - E AtjBjXj} |. (11) 
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Using \P t xj | < c y -|y| for ;' E J, as well as (11) and the causality of the 
Aij and the Bj, we have 

|P T x/| <c(|y| + |E PrAuBjXjl) 
\ \j(=J 1/ 

<c(|y|+ L {PAljPrBjPrXA) 

<c(b|+75 E |Pr^|) 
V J<EJ I 

<c(l + 7«I cy) \y\. 

Let ui,($>2, ■ • • ,w n be defined by «i = c and 

co, = c ( 1 + 7^ L wyj, i = 2,3, . . . ,n. 

We have shown that given x, t, and i, |P T x,| < d,-|y| for some d t E 
{wi,W2» • • • > w nl- Since a>„ = max,co,, we have 

t |P T x,| 2 <nu> n 2 |y| 2 forallr>0, 

which shows that x E L n and that |x | is suitably bounded in terms of 
|y | . This completes the proof for the case in which = 0. 

The proof for the > case is similar. Using primarily (7) and the 
hypothesis that A E S^, we find that 

0\P T B k x k \* < \P T B k x k \ ■ \P T y k \ 
for some k. Therefore, 

\P T B k x k \ Zfi-i\y\ (12) 

for some k. By proceeding essentially as indicated above, we can show 
that 

\P T BiXi\<Sl n \y\ for all i and r > 0, (13) 

in which ft n is the number defined by Ui = 0~ l and 



for i = 2,3, . . . ,n. 
From (6) and (13), 



\P T Xi\ <\P T yi\+Z\AijPrBjXi\ 

j 

< |y|+E70»|y| 

j 

< (1 + nytl„)\y\ 
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for each i and all t > 0. Since this shows that x £ L and that \x | is 
bounded as required, our proof is complete. 

3.1 Comments 

To see that U»;>o St, contains the familiar set of nonnegative operators 
defined on a space of vector-valued functions, let H be any mapping of 
L n into itself such that 

t(w it (Hw)i)><r\w\ 2 , wEL", (13) 

in which a is a nonnegative constant. From (13) it is clear that for each 
w £ L n with |u;| ^ there is an index k such that \wk\ 9* and 
(wk, (Hw)k ) > fi\w\ 2 in which = an~ l . Since | w \ 2 > \ Wk | 2 , we observe 
that H^Sp. 

Theorem 1 is of course a result concerning the L-boundedness of so- 
lutions of (4).* By modifying the hypotheses and proof of Theorem 1 
in a direct manner, an analogous result can be obtained concerning the 
L -continuity of solutions (i.e., concerning the L-boundedness of the 
difference (x a — Xb) of a solution x a of (4) that corresponds to y = y a and 
a solution Xf, that corresponds to y = yb with (y a — y*,) £ L). With regard 
to the necessary modifications of the hypotheses concerning A, the fol- 
lowing definition, in which /? > 0, plays a central role. 
Definition of Tp\ Tp is the set of operators H in L n with the following 
property: For each u and u in L n such that | u — v \ ^ 0, there is an index 
k such that \uk — Vk\ ^ and (uk — Vk, (Hu)k — (Hu)k) ^ (3\uk — 

v k \ 2 - 

In order to be more explicit, let (A. 1') denote the assumption that x a 
+ ABx a = y a and Xb + ABxb = yb in which each Ay and Bj are causal 
operators in K and x a , x^ y a , and yb belong to K. Let A.2' be the hy- 
pothesis obtained from A.2 by replacing "(Bjiv,w) > a\w\ 2 for w £ L 
and all j" with "(23/1* — BjV,u — u) > a\ u — v \ 2 for u and v in L and all 
/ ," and let A.3' and H.l ' be the hypotheses obtained in a similar manner 
from A.3 and H.l, respectively. 

Our L -continuity result (whose proof is omitted) is the following. 
Theorem 2: Let H.l' and A.l' through A.3' be satisfied, let (y a — y b ) £ 
L, and let A £ T$ with a + /3 > 0. Then (x a — Xb) £ L, and there is a 
nonnegative continuous function p defined on [0, °° ) that depends only 
on A and B such that p(0) = and \x a — x^l < p(\y a — yb\)- 



* Results along the lines of Theorem 1 for cases in which B is more general than assumed 
here but both A and B are nonnegative operators are given in Ref. 8, where the stability 
of interconnected systems in the sense of Section 2.2 is not explicitly discussed. A non- 
negative-operator approach to the stability of interconnected systems, as well as its relation 
to other approaches, is discussed in Ref. 7. 

3036 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1978 



3.2 A Corollary to Theorem 1 

We shall refer to the following two hypotheses. 
H.2: Each B( is a continuous mapping in L that maps the zero element 
into itself, and there are positive constants c x and c 2 such that for 
alii 

(B £ u,u)> Cl |B £ u| 2 (14) 

\BiU-Biv\ <c 2 \u-v\ (15) 

for u and v in L. 

H.3: K = E 2t L = L 2 , and for each t > P T is the operator associated with 

.E 2 in the example given in Section 2.1. 

There are many cases in which (14) holds* for all u for some c\ > 0, 
but there is no positive a such that (2) is satisfied for all w.+ On the other 
hand, it is clear that there is a positive C\ with the property that (14) is 
met when (2) holds with a > for all w and there is a positive constant 
8 such that |B,w| < 8\w\ for all w. 

Corollary 1: Let H.2 and H.3 (as well as A.l and A.3) be satisfied. If A 
E S , then (4) is L-stable. 

Proof: Assume that the hypotheses of the corollary are satisfied and let 
/ and /„, respectively, denote the identity operators in K and K n . With 
regard to the following lemma, two elements u and v ofE 2 are taken to 
be the same if and only if |P T (u - v )| = for all r > 0. 
Lemma 1: Let H.3 hold, let Fbea continuous mapping ofL 2 into itself 
such that for some positive constant c < 1 we have 

\Fu — Fv\ < c\u — v\ for u and v in L 2 , (16) 

and let F also be a causal mapping ofE 2 into E 2 . Then (I - F) _1 exists 
and is causal on both L 2 and E 2 . 

Proof of Lemma 1 : Let the hypotheses of the lemma be met. In view of 
(16) and the continuity of F, the equation x - Fx = h with h E L 2 has 
in L 2 a unique solution x which is given by x = lim„^a, x (n) in which x (n) 
= h + Fx n ~ l for n > 1 and x< 0) = h. Thus, (/ - F)~ l exists on L 2 , and 
since h + F(«) is causal on L 2 so is (/ — F)~ l . 

Now let h E E 2 , and for each r > let z T be the unique element of L 2 
that satisfies z T - Fz T = P T h. Since (/ - F)" 1 is causal on L 2 , it is clear 
that P n z r 2 = PriZn for t 2 > ri. Let x be the element of E 2 defined by the 
condition that P^ = P T 2 r for all t > 0. For any t > 0, P^c - P T Fx = P T ? T 
- P T FPrX = P T z T - P t FPtZ t = P t Zt- P t Fz t = PA Therefore x satisfies 
x -Fx = h. Suppose that x x and x 2 in E 2 satisfy h = xi - Fxi = x 2 - Fx 2 



* This type of inequality is among those used in Ref. 8. 

♦ We mention two simple examples: Let L = L 2 and let B, be defined by the condition 
that for each t > 0, (Bnv)(t) = w{t) for \w(t)\ < 1 and (Biiv)(t) = sgn(w(t)) for \w(t)\ > 
1. Then (14) with ci = 1 holds for all u £ L 2 , but there is no a > for which (2) is satisfied 
for all w E ^2- It is not difficult to show that a similar conclusion is reached when B, is the 
convolution operator in L 2 with impulse response e~ l . 
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with \P T (xi — *2)| 5^ for some t > 0. Since F is causal, we have P T hi 
= P T fi2 in which h\ = P T x\ — FP T x\ and h<i = P T X2 — FPjX^ This con- 
tradicts the fact that (/ — F)~ l is causal on L2. Thus, x is the unique 
solution in E2 of x — Fx = h, which means that (/ — F) -1 exists on E2, 
because h is arbitrary. In view of the fact that the solution x of x — Fx 
= h satisfies P T x = P T z T where z T — Fz r — P T h for every t > 0, it is evi- 
dent that the operator (/ — F) _1 on E2 is causal. This proves the 
lemma. 

Let c be a positive constant such that c < min (ci,C2 _1 ). By Lemma 
1, (J„ — cB)~ l exists on K n and B{I n — cB)~ l is causal on K n and maps 
L n into itself. In particular, the equation x + ABx = y can be written 
as 

h + (A + d n )B(I n - cB)-^ = y 
in which h = x — cBx. From A G S , it follows at once that (A + cI n ) 

es c . 

Also, from h = x — cBx and the fact that B is causal on K n and satisfies 
|J5u| <c 2 |u| for a E £", with cc 2 <1, we have IP^ I < (1 -cc 2 )~ 1 |P T A| 
for t > 0. 

Therefore, by Theorem 1, to complete the proof it suffices to observe 
that for any w £: L n , 

(B(I n - cB)~ l w,w) = (Bu,u - cBu) 

= (Bu,u) ~ c\Bu\ 2 

>0 

in which of course u = (I n — cB)~ l w. 

IV. RESULTS CONCERNING THE MATRIX CASE 

Of importance in the theory of interconnected systems is the special 
case in which A is represented by a real nXn matrix a with elements ay, 
in the sense that for each i, 

n 

{Aw)i = £ aijWj, w E £"• 

Throughout this section, "A G M" means that A has such a represen- 
tation with representation matrix a, and, assuming that A^M, Uq(U) 
denotes the set of representation matrices such that A G So (A £ Sp 
with /? > 0). In addition, Po(P) denotes the set of real square matrices 
with nonnegative (positive) principal minors. 

Proposition 1: If (A.l through A.3 are satisfied and) A 6M with a (£ 
Pq, then (4) is not L-stable for some B. 

Proof: Let the hypotheses be met, and let 1„ denote the identity matrix 
of order n. From a (£ Po, it follows that there is a diagonal matrix d = 
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diag (di,d 2 , • ■ . ,d„) with d, > for all i such that (d + a) and hence 
(1„ + ad~ l ) are singular.* 

For each i, let B, be defined by the condition that BiW = d, " l w for w 
G K. Let u be any real nonzero n -vector that is annihilated by (1„ + 
ad -1 ), and let e be any element of K different from the zero element 
of L. With x = ue, we have x + ABx = d. This shows that (4) is not In- 
stable for the particular B constructed. 

In order to proceed to the first theorem in this section, it is necessary 
to introduce the following definitions. For each w E L n , a w denotes the 
matrix obtained from a by replacing a, ; with a(j(Wi,Wj) for all i and /. For 
any w E L n with \w\ 5^0, T{a w ) denotes the matrix obtained from a w 
by deleting both the ith row and ith column for all i E \i' \wt\ = 0}. 
Theorem 3: Let A^M.We have a E U(a E U ) if and only if T(a w ) 
E P(a w E Po) for each w E L n with \w\ ^ 0. 

Proof of Theorem 3: We shall use two lemmas. With regard to the first 
of the lemmas, M n denotes the normed linear space of real nXn matrices, 
with the usual Euclidean norm, and C denotes |uGM n : there is a w E 
L n with \w\ = 1 such that u (J = (w it Wj) for all i and;'}. 
Lemma 2: C is compact. 

Proof of Lemma 2: The set C is obviously bounded. To show that C is 
closed, let u (1) ,u (2) , ... be a sequence of elements of C that converges to 
some element u of M n . 

Given a real n -vector u, for each w E £" and its corresponding element 
u of C, we have v tr uv = (2fDft0{,2j!>£i^) > O.t Thus, each i*W is nonne- 
gative definite, and therefore it follows that u is nonnegative definite. 
In view of the fact that Tr(u { J ] ) = 1 for all;', it also follows that Tr{d) = 
1. 

Since u is nonnegative definite and has unit trace, there is an or- 
thogonal matrix T and a diagonal matrix D = diag (di,d 2 , ■ ■ ■ ,d n ) with 
dj > and 2yd/ = 1 such that 

t 2 = TDT"" = Ld/T/(T,)"-, 

/ 

in which T t is the /th column of T. Referring to the pairwise mutually 
orthogonal elements ei,e 2 , . . . ,e n mentioned in Section 2.1, let z in L n 
be defined by 

z = £d ( i/27> ( -. 
t 

Using the orthonormality of the e,, and the fact that for each / the sum 

of the squares of the components of T\ is unity, it is not difficult to verify 

that \z\ = 1, and that (Z|,zy) is equal to the ij'th element of u for all i and 



* A proof is given in Ref. 10. Another proof can be obtained from the fact that, since a (£ 
Pa, there is (see Ref. 9) a real nonzero n -vector q such that q, (aq)i < for every i such that 
qt * 0. 

+ The superscript "tr" denotes transpose. 
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j. This shows that C is closed, and completes the proof of the 
lemma.* 

The following lemma is proved in Ref. 9. 
Lemma 3: A real square matrix m belongs to P(Po) if and only ifuk(mu)k 
> (vk j* and Vk (mu)k ^ 0) for some k for each real nonzero vector 
v of dimension equal to the order of m. 

In order to prove the theorem, suppose initially that T(a w ) £ P for 
every w £ L n with |w| ^ 0. By Lemma 3, for each w £ L n with \w\ = 
1 we have (T(a w )v)k > for some index k, when all of the components 
of the vector v of compatible dimension are unity. Thus, max; 2" =1 
aij(wi,Wj), which we view as a function of the matrix u whose elements 
are the (Wi,wj), is positive for each w in L n with unit norm. Since max t 
2"= i aij(wi,Wj) is obviously a continuous function of u, and, by Lemma 
2, C is compact, there is a a > such that 

n 

minmax £ aij{Wi,Wj) = a. (17) 

Therefore, for each w £ L n with |u;| = 1 there is an index k such that 

n 

£ a kj (w k ,Wj) > a\w\ 2 > a\w k \ 2 , 

from which we see that for each w €EL n with \w\ ^ 0, there is a k such 
that | u>k | t* and 

E ajv(l0A,IO;)£*|u>ft| S . 

Thus, a e U. 

To show that oGt/o when a„, G Po (and hence r(au,) £ Po) for each 
w £ L n with |«; | > 0, we observe that then, by Lemma 3, for each w £ 
L n with \w\ 9^ we have (r(au,)v)A > for some k when the components 
of v are all unity. Therefore, for each w £ L" with | u; | 5^ 0, there is a fe 
such that \u>k\ 7* and 

n 

L a kj (w k ,Wj) > 0, 

which means that a £ (/o- 

Suppose now that for some w £ L n with |u>| > we have T(a w ) (£ 
P(Ou; € Po)- Then, by Lemma 3, there is a nonzero vector v such that 
Uk(T{a w )u) k ^ (uk(T(a w )v)k < 0) for every k such that u* 5^ 0. Thus, 
by multiplying each wi for which \u>i\ ^ by the appropriate component 
of v, it is a simple matter to construct a 2 £ L n for which |z | 5^ and 
Xjaij(zi,Zj) < for all ftZ/tyfai*/) < for all 1 such that |z,- 1 ^ 0). This 
completes the proof of the theorem. 

Corollary 2: Ifn>3, U( Uq) is a proper subset of the matrices of order 
n in P(P ). 



* Of some peripheral interest is the fact that it is not necessary to assume that L is 
complete. 
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Proof: To see that U(U ) is a subset of P(P ), let a £ l/(L/ ), let e be any 
element of L such that | e | = 1, and let w be the element of L n defined 
by wi = e for all i. Thus a w = a, and, by Theorem 3, a £ P(Po)- 

In order to show that for n > 3 there is a matrix of order n in P(Po) that 
is not contained in U(Uo), observe that it is sufficient to consider the 
n = 3 case, and let a (+) and a (0) be defined by 



( (+) = 



,(0) _ 



1.1 1 - 


10 


1 1.1 1 


1 1 1.1, 


1 1 -101 




1 1 1 




1 1 1. 





We have a<+> £ P and a<°> £ P . Let w £ L 3 be given by 




in which e\ and e<i are orthogonal elements of L with unit norm. It is a 
simple matter to verify that 





[1.1 


l 


-10 


a (+) = 


1 


2.2 


11 




, 1 


11 


111.1 




ri 


1 


-10] 


a (0) = 


l 


2 


11 






l 


11 


101 



and that we have det[a, ( + ) ] < and det[a j„ 0) ] < 0. By Theorem 3, this shows 
that a (+) (a (0) ) 6 U(U ), which completes the proof. 
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Corollary 3: If n = 2, U(U ) = P(P ) restricted to 2X2 matrices. 
This follows directly from Theorem 3.* 

4.1 Definitions 

Let N denote the set of real symmetric nonnegative definite matrices 
m of order n with ma > for some i. For each mEW, let a m denote the 
nXn matrix whose ijth element is a^mij for all i and /,+ and let A(a m ) 
denote the matrix obtained from a m by deleting row i and column i for 
each i £ [/": rrijj = 0}. 

Corollary 4: We have a £ L/(l/ ) if and only if A(a m ) £ P(a m £ P ) for 
each m £ N. 

Proof: The proof of Lemma 2* shows that a real matrix m of order n 
belongs to N if and only if there is a w in L" such that | w \ ?£ and my 
= (wi.Wj) for all i and y. Thus, Corollary 4 follows from Theorem 3. 

4.2 Introduction to Corollary 5 

In order to present our next corollary, we need the following additional 
definitions: Let S(m) denote the set of all matrices obtainable from a 
given real nXn matrix m by replacing each off-diagonal element m^ of 
m with rijiriij, where the r# are real numbers that satisfy ry = r, t and fr# | 
< 1. Let R denote jm £ P: S(m) C P}, and, similarly, let R = [m £ P : 
S(m) C Po}. 

When m = 2 and P and Po are restricted to 2X2 matrices, we have R 
= P and R = P . On the other hand, if we let a< + >(\) and a<°>(\), re- 
spectively, denote the matrices obtained from a (+) and a (0) of the proof 
of Corollary 2 by multiplying the (1,3) and (3,1) elements by a scalar 
variable X, then a( + >(l) £ P and a«»(l) £ P , but a<+>(0) € P, and, 
similarly, a (0) (0) (£ P . This shows that R(Ro) is a proper subset of the 
nXn matrices in P(Po) when n > 3. s Two familiar classes of matrices 
contained, for example, in Rq are the set of row-sum dominant matrices 
and the set of column-sum dominant matrices. 
Corollary 5: We have a £ U(Uo) if either 
(i) a £ R(R ). 

(ii) There are diagonal matrices d\ and d% of order n with positive di- 
agonal elements such that d±ad2 is positive definite {nonnegative 
definite).^ 

Proof: Suppose first that a £ R{Rq). Let w be any element of L n such 
that |«;| 5^ 0, let c = diag(c 1( C2, . . . ,c n ) in which for all i, c t = if \w(\ 



* The proof of Corollary 2 shows that U C P and U C Po for n > 2. 

f In other words, let a m denote the "Schur product" of a and m. 

I Lemma 2 is used in the proof of Theorem 3. 

s It will become clear that this proposition also follows from Corollary 2 and Corol- 
lary 5. 

" As usual, we say that a real square matrix m is positive definite (nonnegative definite) 
if and only if the symmetric part of m is the matrix of a positive definite (nonnegative 
definite) quadratic form. 
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= and c, = \w l \ _1 if \w ( \ > 0, and let r(ca w c) denote the matrix ob- 
tained from (ca w c) by deleting the rows and columns corresponding to 
the indices i for which c, = 0. In view of the fact that | (w it Wj) \ < \wj \-\u)j\ 
for each i and;', we see that T(ca w c) E P(Po) and hence* that T(a w ) E 
P(P ). By Theorem 3, a E U(U ). 

At this point, we need the following lemma. 
Lemma 4^: Ifp = {pij\ and q = \qij\ are real square matrices of the same 
order, with p positive definite and q symmetric, nonnegatiue definite, 
and such that qu > for all i, then r = {p,,*?,;! is positive definite. 
Proof of Lemma 4: Let p and q be as indicated, and let k denote the order 
of p. The proof of Lemma 2 shows that L 2 contains k functions 
/11/21 ■ ■ ■ >fk such that 

Qij= f~fi(t)fj(t) dt for alii and;'. 

With v any real nonzero k -vector and with X the smallest eigenvalue of 
the symmetric part of p, we have 

v tr rv = Y.ViVjPij f fi(t)fj(t)dt 

ij •'O 

= ("HPijUifMVjfj^dt 
JO i.j 

^ r ro xE(^(o) 2 ^ 

Jo , 

>0, 

which shows that r is positive definite. 

To complete the proof of the corollary, suppose that d\ad2 is positive 
definite, with d x and d 2 as described, and let m (E N. By Lemma 4, 
A\dia m d 2 ] (i.e., A{a m ) with a replaced with diad 2 ) is positive definite and 
hence it belongs to P. Therefore, A|a m ) G P, and, by Corollary 4, 

aE u. 

The proof for the case in which d\ad 2 is nonnegative definite is es- 
sentially the same, and is omitted. 

4.3 Comments Regarding Corollary 5 

In light of the fact that R = P restricted to 2X2 matrices when n = 
2, the following special result is a direct consequence of Corollaries 5 and 
1, and the content of the proof of Proposition 1. 
Proposition 2: Letn = 2 and A E M. Let H.2 and H.3 {as well as A.l and 



* Here and in another part of the proof, we use the easily proved result that a real square 
matrix m belongs to P(P ) if and only i(dimd 2 £ P(Po) for every pair of compatible di- 
agonal matrices d x and d 2 with positive diagonal elements. 

f A proof that the conclusion of Lemma 4 holds when q is positive definite is given in 
Ref. 11. 
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A. 3) be satisfied. Then (4) is L-stable for every B if and only if a £ 

An example of a matrix a that is nonnegative definite and such that 
a (£ Rq is given by 



a = 



1 1 1 
1 2 1/2 2 1/2 
1 2 1 ' 2 2 1/2 



since 



det 



1 


1 


1 


1 


2 l/2 


2 l/2 





21/2 


2 l/2_ 



<0. 



Similarly, a very simple example of an a £ Ro such that diad 2 is non- 
negative definite for no suitable di and d 2 is 



Lo oj 



Of some interest is the fact that a £ i/ if a is an M-matrix (i.e., if a 
has positive principal minors and nonpositive off-diagonal elements); 
in that case there is a diagonal matrix d with positive diagonal elements 
such that ad is strongly row-sum dominant* and therefore ad and con- 
sequently a belong to R. 
Theorem 4: Let A (EM. If a,, = for all i, then a £ Uq if and only if a 

GRo- 

Proof: The "if part" is a special case of Corollary 5. 

Suppose that a £ £/n with a„ = for all i, and suppose also that a (£ 
Ro in which case there is an element b of S(a) such that b € Po- Let b 
be given by 6,7 = for all i, and 6 i; = r^aij with r, ; ■ r# for i 9^ j. Choose 
n real numbers rn,r 22, ■ ■ ■ Jnn so that the nXn matrix m given by ttty 
= Ty for all i and / is nonnegative definite. Observe that m £ N. Since 
a m = fe, by Corollary 4, we have a contradiction to the supposition that 
£ £/o- Therefore, a £ /?o when a £ t/o and a,, = for all i, which 
completes the proof of the theorem. 

4.4 Comment regarding Theorem 4 

We can have a £ Po with a„ = for all i, and a € R . For example, 
let 



a = 



[0 


-1 1] 


1 





.0 


1 0. 



* See the theorem given on page 387 of Ref. 12. 
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Fig. 1 — Flow graph of an interconnected system. 

Then a E Po, and 



det 



rO -r 12 r 13 
ri 2 
r 23 J 



- ^23^*12^13 < 



for, say, r 23 = r l2 = -r i3 = 1. 

4.5 A Specific example of an L-stable Interconnected system 

Assume that H.2 and H.3 (as well as A.l and A.3) are satisfied. 
For the system described in flow-graph form in Fig. 1, we have 

y\ = *i 

y 2 = x 2 - BiXi + B3X3 - B4X4 
y 3 = x 3 - 2B 2 *2 + 2B4X4 
y 4 = x 4 + B 2 x 2 - B 3 x 3 . 
Here A G M , with 


-1 

2 



a = 



' 








-1 





1 





-2 





. 


1 


-1 



To see that a £ Rq, consider the matrix 




-r 21 

-2r 23 
, r 24 



^"23 



"^34 





-r 2 4 

2r 34 
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We observe that its determinant vanishes and every IX 1 and 2X2 
principal minor is nonnegative for all real values of the r, ; . It is a simple 
matter to verify that its principal minor of order three obtained by de- 
leting the first row and first column vanishes for all values of the r#, and 
it is clear that every other principal minor of order three also vanishes 
for all values of the ry. 

Since a £ Rq, by Corollary 1 and either Corollary 5 or Theorem 4, the 
system described in Fig. 1 is L -stable. 

Another way to prove that the system in Fig. 1 is L -stable is as follows. 
Since H.2 holds, \B\u\ < C2I u | for u £ L. It therefore suffices to show 
the L -stability of the system obtained from the flow graph in Fig. 1 by 
deleting B\, x\, andyi. That can be done with the aid of Corollary 5 by 
verifying that the interconnection matrix a of the modified system has 
the property that there is a 3X3 diagonal matrix d with positive diagonal 
elements such that da is nonnegative definite. 
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